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ELECTRON-POSITRON-PHOTON PLASMA AROUND A 
COLLAPSING STAR 



R. RUFFINI, L. VITAGLIANO AND S.-S. XUE 

We describe electron-positron pairs creation around an electrically charged star 
core collapsing to an electromagnetic black hole (EMBH), as well as pairs annihi- 
lation into photons. We use the kinetic Vlasov equation formalism for the pairs 
and photons and show that a regime of plasma oscillations is established around 
the core. As a byproduct of our analysis we can provide an estimate for the ther- 
malization time scale. 



1. Dynamics of Dyadosphere 

Dyadosphere was first introduced in Ref. 1 as the region surrounding an 
electromagnetic black hole (EMBH) in which the electromagnetic field 
strength exceeds the critical value Sc for electron-positron pair creation via 
the mechanism a la Heisenberg-Euler-Schwinger.^''^ The relevance of the 
dyadosphere around an EMBH, for the astrophysics of gamma-ray busts 
has been discussed in Refs. 1, 4-6 (the external radius of dyadosphere will 
be denoted by r^s)- In those papers the pair production in dyadosphere has 
been described as an electrostatic problem: instantaneously a massive body 
collapses to an EMBH whose charge is large enough that the electric field 
strength £ exceeds £c and the Schwinger process is triggered in the entire 
dyadosphere; moreover the pairs are produced at rest and remain at rest 
during the whole history of their production; finally they instantaneously 
thermalize to a plasma configuration (see Fig. 1). These ansatz, formu- 
lated for the sake of simplicity, allow one to estimate the number density 
of pairs produced as well as the energy density deposited on the pairs in 
a straightforward manner. We relax the hypothesis that the large electric 
field is instantaneously built up and take the following dynamical point of 
view: 

(1) A spherically symmetric star core endowed with electric, say posi- 
tive, charge Q, collapses. We assume that the electromagnetic field 
strength £ on the surface of the core is amplified to £c during the 
collapse and the Schwinger process begins. 
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Figure 1. Vacuum polarization process of energy extra<;tion from an EMBH. Pairs axe 

created by vacuum polarization in the dyadosphere and the system thermalizes to a 
neutral plama configuration (see Ref. [1] for details). 

(2) The pairs produced by the vacuum polarization progressively screen 
the electromagnetic field of the core, thus reducing its strength. 
Furthermore the charges (electrons and positrons) are accelerated 
by the Lorentz force in the electromagnetic field. Finally particles 
and antiparticles annihilate into photons. 

An enormous amount of pairs {N ~ ^ aH'' ^'^^'^^ '^C is the Compton 
length of the electron) is produced, as claimed in Rcfs. 1, 4-6, if the core 
charge is not annihilated by the charge of the accelerated electrons during 
the gravitational collapse (see Ref. 7). Therefore it is useful to study the 
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dynamics of the electron-positron-photon plasma in the electric field of the 
core in some details. This will be the main object of the next section. As a 
byproduct of the analysis we obtain an estimate for the time scale needed 
for the thermalization of the system. 

In Ref. 8 it was suggested that the exact solution of Einstein-Maxwell 
equations describing the gravitational collapse of a thin charged shell can 
be used as a simplyfied analytical model for the gravitational collapse of 
a charged core; it was also discussed in some details the amplification of 
electromagnetic field strength on the surface of the core. Here we briefiy 
review some of the results of Ref. 8. The region of space-time external to 
the core is Reissner-Nordstrom with line element 

in Schwarzschild like coordinate (t, r, 6, </>) , where dfl^ = dO"^ + cos^ Odcp'^ , 
f — f [r) — 1 — ^ + ^ ; M is the total energy of the core as measured at 
infinity and Q is its total charge. Let us label with R and T the radial and 
time-like coordinate of the shell, then the equation of motion of the core is 
(cfr. Ref. 8) 

If = -|^v/F^(SF7CR) (2) 

where H (i?) = ~ ^jMofl i -^o is the rest mass of the shell. The analytical 
solution of Eq. (2) was found in Ref. 8 in the form T = T (R). According 
to a static observer O placed at the event xq = {R, T {R) ,6o,(po) the core 
collapses with speed given by 

where dR* = f~^^'^dR and dT* = f^/'^dT are spatial and temporal proper 
distances as measured by O. In Fig. 2 we plot V* as a fmiction of R for a 
core with M ^Mo = 2OM0^ miA £, = ^ = IQ-^, IQ-^, IQ-^. Recall that 
dyadosphere radius is given by^ 

rds = (4) 

where c is the speed of light; e and m,, are electron charge and mass re- 
spectively. Then note that V^^ = ~ 0.2c for ^ = 0.1. 



"The condition Mq = M corresponds to a shell starting its collapse at rest at infinity. 
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Figure 2. Collapse velocity of a cliargcd star of mass Mq = M = 20Mq as measured 
by static observers as a function of the radial coordinate of the star surface. As the 
charge is not too large < 0.1) there is not much difference between collapse velocities 
of stars with different charge. Dyadosphere radii for different charge to mass ratios 
= 10~^, 10~^, 10~^) are indicated in the plot together with the corresponding velocity. 

2. Plasma oscillations and screening 

Wc now turn to the pair creation taking place during the gravitational 
collapse. The gravitational fields of the core is considered classical; the 
gravitational effects of the electron-positrons-photons plasma are neglected. 

The most detailed framework for studying electromagnetic vacuum po- 
larization and particlc-antiparticlc scattering around an electromagnetic 
collapsing core is quantum electrodynamics in the classical external elec- 
tromagnetic field of the core on the Reissner-Nordstrom space-time around 
the core itself. Of course a number of approximations is needed in order to 
make the problem be tractable. Let us discuss such approximations. 



(Homogeneity) First of all, the static Reissner-Nordstrom space-time region exter- 



nal to the collapsing core is naturally splittcd in space (hypersur- 
faces orthogonal to the static Killing field) and time. In the local 
frames associated with static observers, the electromagnetic field of 
the core is purely electric. Moreover, we will see that the length 
scale L over which the electric field as well as the particle num- 
ber densities vary, is much larger than the length scale I which is 
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charachteristic of the electron-positron motion. Thus we can divide 
dyadosphere into smaU regions Di 

Di:ri<r< n+i = + e; (5) 

r-+ < Ti < Tds £<l 

such that for any i the system formed by the electric field and the 
pairs can be considered homogeneous in Di. 
(Flat space-time) For, in geometric units, the electron charge e is miich larger than 
the electron mass nie, the gravitational acceleration is negligible 
with respect to electric acceleration for sufficiently large electric 
field strengths (c;ven much less than £c), therefore wc will neglect the 
curvature of space-time and use the local frame of a static observer 
as a globally inertial frame of the Minkowski space-time. 
(Mean field) The number of pairs is so high that a semiclassical formalism and 
mean field approach can be used, in which the total electromagnetic 
field (core electromagnetic field and screen field due to pairs) is con- 
sidered to be classical, while the electron-positron field is quantized. 
It has been shown^"^^ that, if we neglect scattering between parti- 
cles, the semiclassical evolution of the homogeneous system in a flat 
space-time is well described by a Boltzmann-Vlasov-Maxwell system 
of partial differential equations, where the electrons and positrons 
arc described by a distribution function /,. = /,. (t, p) in the phase 
space, where t is the inertial time and p the 3— momentum of elec- 
trons. Finally we use the method presented in Ref. 13 to simplify 
such a Boltzmann-Vlasov-Maxwell system. 

Let us summarize results in Ref. 13: we obtained the following system of 
ordinary differential equations which simultaneously describes the creation 
and evolution of electron-positron pairs in a strong electric field as well as 
the annihilation of pairs into photons: 

' ^rie = S(£)- 2nlaipj^ \Tre\\ \ + 2n^(72 
An-y = 4n2(Tip-i [TTeii I - 4n^<T2 

JjPe = erieSpe^ [TTell | 4" ^Sjp {£) - 2neP^CFip~^ \TTe\\ I + 2n^p^(J2 .g-. 
j-^p^ = AUePeO-lPe^ \^e\\ \ - ^n^p^CT2 
^TTell = eUeS - 2ne1Te\\0-lP~^ \TTe\\ \ 

, ft^ = -ZeneP^i iTTell | - jp {£) 

where Ue {n-^ ) is the electron (photon) number-density, p,, (p, ) is the elec- 
tron (photon) energy-density, TTgy is the density of electron radial momen- 
tum and £ the electric field strength. Finally S {£) is the Schwinger prob- 
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ability rate of pair creation, jp {£) is the polarization current-density, ai 2 
are total cross sections for the processes e+e" 77 and the corresponding 
terms describe probability rates of pair annihilation into photons and vice 
versa. System (6) was numerically integrated in Ref. 13. 

Here it has to be integrated once for each of the regions Di (see (5)) 
with initial conditions 

rie = rij = Pe = Pj = TTe\l =0; £0 = ^. (7) 

Let us recall the main results of the numerical integration. The system 
undergoes plasma oscillations: 

(1) the electric field oscillates with lower and lower amplitude around 

0; 

(2) electrons and positrons oscillates back and forth in the radial direc- 
tion with ultrarelativistic velocity; 

(3) the oscillating charges arc trapped in a thin shell whose radial di- 
mension is given by the elongation Al = \l — lo\ of the oscillations, 
where Iq is the radial coordinate of the centre of oscillation and 

I = 1^ --^dt. (8) 

Note that ^ = is the radial mean velocity of charges (we plot 

the elongation Al as a function of time in Fig. 3); 

(4) the lifetime At of the oscillation is of the order of 10^ — IO^tq (see 
Fig. 3). 

(5) in the time At the system thermalizes in the sense that both num- 
ber and energy equipartition between electron-positron pairs and 

photon are approached. 

In Fig. 4 we plot electrons mean velocity v as a function of the elon- 
gation during the first half period of oscillation, which shows precisely the 
oscillatory behaviour. 

3. Conclusions 

In a paper under preparation^ we are examining the conditions under which 
the charge of the collapsing core is not annihilated due to vacuum polar- 
ization as a consequence of the above plasma oscillations. 

Note that e+e~ ^ 77 scatterings is marginal at early times (t <C 
At) since the cross sections ai^2 are negligible in the beginning of pair 
production. However at late times {t > At) the system is expected to 
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Figure 3. Electrons elongation as function of time in the case r = ^f'ds- The oscillations 
are damped in a time of the order of IO'^tq- 

ic;lax to a plasma configuration of thermal equilibrium.^^ Thus a regime of 
thermalized electrons-positrons-photons plasma begins in which the system 
can be described by hydrodynamic equations. It is shown in Refs. 12, 14 
that the equations of hydrodynamic imply the expansion of the system. 
In "brief" the system reaches the ultrarclativistic velocities required in a 
realistic model for GRBs. It is worthy to remark that the time scale of 
hydrodynamic evolution {t ~ 0.1s) is, in any case, much larger than the 
time scale At needed for thermalization. 
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Figure 4. Electrons mean velocity as a function of the elongation during the first half 
oscillation. The plot summarize the oscillatory behaviour; as the electrons move, the 
mean velocity grows up from to the speed of light and then falls down at again. 



8. C. Cherubini, R. Ruffini and L. Vitagliano, Phys. Lett. B545, 226 (2002). 

9. Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper and E. Mottola, Phys. 
Rev. Lett. 67, 2427 (1991). 

10. Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper and E. Mottola, Phys. 
Rev. D45, 4659 (1992). 

11. F. Cooper, J. M. Eisenberg, Y. Kluger, E. Mottola and B. Svetitsky, Phys. 
Rev. D48, 190 (1993). 

12. R. Ruffini, J. D. Salmonson, J. R. Wilson and S.-S. Xue, A&A 359, 855 
(2000). 

13. R. Ruffini, L. Vitagliano and S.-S. Xue, Phys. Lett. B559 (2003) 12. 

14. R. Ruffini, L. Vitagliano and S.-S. Xue, (2003) in these proceedings. 



